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ABSTRACT 

Deep Very Large Array imaging of the quasar 3C 345 at 4.86 and 8.44 GHz has been used to 
study the structure and linear polarization of its radio jet on scales ranging from 2 to 30 kpc. There 
is a 7-8 Jy unresolved core with spectral index a ~ —0.24 [I v oc v a ). The jet (typical intensity 
15 mJy/beam) consists of a 2.5" straight section containing two knots, and two additional non-co- 
linear knots at the end. The jet's total projected length is about 27 kpc. The spectral index of the 
jet varies over —1.1 < a < —0.5. The jet diverges with a semi-opening angle of about 9°, and is 
nearly constant in integrated brightness over its length. A faint feature north-east of the core does 
not appear to be a true counter-jet, but rather an extended lobe of this FR-II radio source seen in 
projection. The absence of a counter-jet is sufficient to place modest constraints on the speed of the 
jet on these scales, requiring /3 > 0.5. Despite the indication of jet precession in the total intensity 
structure, the polarization images suggest instead a jet re-directed at least twice by collisions with 
the external medium. 

Surprisingly, the electric vector position angles in the main body of the jet are neither longitudinal 
nor transverse, but make an angle of about 55° with the jet axis in the middle while along the edges 
the vectors are transverse, suggesting a helical magnetic field. There is no significant Faraday rotation 
in the source, so that is not the cause of the twist. The fractional polarization in the jet averages 
25% and is higher at the edges. In a companion paper dRoberts & Wardle||2012|) it is shown that 
differential Doppler boosting in a diverging relativisitic velocity field can explain the electric vector 
pattern in the jet. 

Subject headings: galaxies: active — galaxies: jets — quasars: individual (3C 345) 



1. INTRODUCTION 

The compact radio source 3C345 (J1642+3948) is a 
16th magnitude quasar at a redshift of z = 0.593 (6.5 
kpc/arcsecond in concordance cosmology, h = 0.73). 
This is one of the best studied quasars at a wide range of 
wavelengths and angular resolutions because it is nearby 
and bright F] The source is variable in X-ray, optical, 
and radio bands on time scales ranging from minutes 
to months, and i t is highly polarized at optical and r a- 
dio wavelengths QKinman et al| ( |1968D, |Aller| 11970^). 
The radio structu re was studied by |Kollgaard, Wardle, 
|fc Roberts] (11989) on angular scales from 0.3 through 
30" (2 to 20CTIpcj using multi-array VLA data at 5 GHz. 
They were able to trace the radio source from the core 
north-west along the jet, which then turns northward and 
eventual ly loops around to the east and then t o the south 
(see also Murphy, Bro wne, fc Perley] ( |1993|) ). |Kollgaard^| 
Wardle, ki Roberts (1989) concluded that this highly dis- 



torted structure was the result of projection, and that 
3C345 was a Fanar off- Riley Type II radio source (Fa- 



naroff fc Riley 1974[ ) seen nearly end-on. 3C 345 has been 
observed extensively by Very Long Baseline Inter ferome- 
try, and is sup erluminal with in ferred Lorentz factor V up 
to at least 20 ( |Lister et al.|2009[ ) . Thus it is assumed that 
the jet near the core is aligned within an angle i ~ 1/Y 
of the line of sight, consistent with the interpretation of 

roberts@brandeis.edu 

1 There are over 1000 references to 3C 345 in the NASA Extra- 
galactic Database (http:/ / nedwww.ipac.caltech.edu/ ), to which we 
refer the reader for further information. 



the large-scale structure. Since at larger distances the 
jet is strongly curving, the inclination of the jet within 
about 0-4" of the core — the region studied below — is not 
well determined. This uncertainty makes it impossible 
to de-project the jet with complete confidence; however, 
it seems certain that bends and opening angles are exag- 
gerated by projection. 

In this paper we use high-resolution, high dynamic 
range VLA images at 5 and 8 GHz to study the structure 
and linear polarization of 3C 345. The observations and 
data reduction are described in ^2| the properties of the 
jet described in ^3j we discuss the results in Q and our 
conclusions are presented in ^5] 

2. OBSERVATIONS AND DATA REDUCTION 

Interferometer data from 1989 January 7-8 (epoch 
1989.02) were obtained from the VLA data archive (ex- 
periment code AR196). The array was in the A con- 
figuration with 27 working antennas. The frequencies 
used were 4.86 GHz (C band, A6.2 cm) and 8.44 GHz 
(X band, A3. 6 cm), with a bandwidth of 50 MHz for 
each of the two IF systems at each band. T he data were 
edited and calibrated in AIPS flNRAO|20lT| ). The nearby 
source 1635+381 was used for phase calibration, and the 
primary flux calibrator 3C 286 was used for amplitude 
calibration. To obtain the maximum dynamic range in 
the images a baseline-dependent calibration was executed 
with the task BLCAL. The baseline calibration tables 
were determined from the data for 3C 286 and its AIPS 
CLEAN component models. The tables were then at- 
tached to the data files of 3C 345 for their respective fre- 
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quency bands. It was found that making a baseline cor- 
rection using an image of the bright source 3C 84 made 
from data collected about 12 hours earlier produced es- 
sentially the same results at 5 GHz, but was superior 
at 8 GHz. Hence this correction was used at the higher 
frequency. Baseline-corrected images had signal-to-noise 
about a factor of two better than those made without 
baseline correction, and were judged to have higher im- 
age fidelity. 

Polarization calibration was done with the phase cali- 
brator using the task PCAL. The position angle of linear 
polarization on the sky was set by 3C 286. We estimate 
that the absolute position angle calibration is good to 
within ±3° at 5 GHz and ±5° at 8 GHz. 

Images were made using AIPS a nd the Caltech dif- 
ference imaging package DIFMAP ( |Shepherd, Pearson^] 
fc Taylor|[l994| as follows. Using DIFMAP, at each 
step in theSELFCAL-CLEAN cycle, a tight CLEAN 
box was placed around the brightest remaining flux and 
purely positive CLEAN components were removed until 
the process ended. Then a 15 minute phase-only SELF- 
CAL was done, and the boxing-CLEANing-SELFCAL 
resumed. The u-v weighting was changed and u-v tapers 
applied as required to recover more flux. Only when this 
process terminated was an amplitude SELFCAL used, 
with a 60 minute timescale. Phase-only SELFCALs were 
resumed until the process terminated. Finally, the self- 
calibrated visibilities were moved into AIPS and the final 
images made with the task IMAGR using the restoring 
beams listed in Table [I] To create images that could be 
compared at the two frequencies, uniform weighting (RO- 
BUST = — 5) was used for 5 GHz and natural weighting 
(ROBUST = 5) for 8 GHz, a u-v taper with parameters 
(850, 800) kA was used at both frequencies, and identical 
circular restoring beams of full-width at half-maximum 
0.30" were applied. To study the divergence of the jet, 
full-resolution images were made from the data at 8 GHz; 
it was judged that "optimum weighting" (ROBUST pa- 
rameter 0) and a CLEAN beam of 0.20" gave the best 
compromise between resolution and signal-to- noise. To 
simplify the analysis and image display, most of the im- 
ages were rotated by —50° in position angle. 

In Table [l] we give the parameters of the images and 
the figures. The root-mean-square (RMS) noise reported 
is from a 5" x 5" box centered 7.5" southeast of the core, 
and is reasonably representative of the uncertainties in 
the images (see below). 

3. PROPERTIES OF THE 3C345 JET AT 5 AND 8 GHZ 
3.1. Total Intensity Structure 

A 20"-wide field of 3C 345 at 5 GHz at its correct orien- 
tation on the sky is shown in Figure [I] The detailed total 
intensity structures at 5 and 8 GHz are shown in Figure [2j 
At the epoch of observation (1989.02), the core of 3C 345 
was essentially unresolved, and had flux densities of 8.1 
and 7.1 Jy at 5 and 8 GHz, respectively. Figure [3] shows 
a full-resolution image of 3C 345 at 8 GHz. 

Since the jet is weak compared to the core, we ex- 
amined the reliability of the images by comparing those 
made independently from the same data by two of the 
authors (DHR & VVM). The RMS difference between 
the pairs of images was about 0.5 mJy/beam at each fre- 
quency, consistent the RMS noise in the images. Ridge- 
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Fig. 1. — Total intensity image of the quasar 3C 345 made from 
5 GHz VLA A-array data, natural weighting (ROBUST = 5), full 
field showing the orientation of 3C 345 on the sky. This image is 
the mean of images made independently by two of the authors from 
the same data. See Table [l] for image data and figure details. 
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Fig. 2. — Total intensity images of the quasar 3C 345 made from 
(a) 5 &; (b) 8 GHz VLA A-array data. The images are the means of 
those made by two of the authors, and have been rotated by —50° 
in position angle. See Table ^ for image data and figure details. 

line intensities at 5 and 8 GHz are compared in Fig- 
ure [4j and are consistent with an RMS difference of 
~ 0.5 mJy/beam. This gives an idea of the dependence of 
the intensity levels in this faint jet on the personal equa- 
tion in the hybrid imaging technique. This of course does 
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Fig. 3. — VLA image of 3C 345 at 8 GHz made with "optimum" 
weighting (ROBUST = 0). The image is the mean of those made 
by two of the authors, and has been rotated by —50° in position 
angle. See Table ^ for image data and figure details. 

not mean that the absolute intensity levels are reliable 
at that level because any calibration uncertainties cancel 
in this comparison. 

3.2. Total Intensity Properties of the Jet 

The jet consists of a westward straight section^] of ap- 
proximate projected length 2.25" (16 kpc), with two 
"knots" that we will call "Kl" and "K2," and two ad- 
ditional bright spots at the west end of the jet that do 
not lie on the same east- west line. To avoid confusion 
with Kl and K2, we will refer to the first of these spots 
as the "knob," denoted "KB," and the final one as the 
"hotspot," or "HS;" see Figure^. Examination of the 
polarization images below suggests that we further sub- 
divide the knob into a northern piece we will refer to as 
"KBN" and a southern piece "KBS." 

3.2.1. Spectral Index 

Figure [5] shows a spectral index image in color over 
total intensity contours, made from from the images in 
Figure [2] The core has a rather flat spectrum with 
a = —0.24, while the spectral index of the jet along the 
ridge-line ranges over —1.1 < a < —0.5 [I v oc v a ). The 
spectrum is significantly steeper along the south edge of 
the jet yet flatter on the north side, and flatter between 
K2 and KB, and between KB and HS. The spectral index 
is plotted along with the ridge-line intensities in Figure [6j 
The error bars for a were determined from the estimate 
of the systematic uncertainties in the intensities derived 
by comparison of pairs of independent images made from 
the same data. 

3.2.2. Is There a Counter- jet in 3C345? 

Kollgaard, War die, fc Robert s| ( |1989 ) detected a pos- 
sible counter-jet with (true) position angle ~ 55°. In the 
naturally- weighted image (Figure [T]) we see the same fea- 
ture in essentially the same location at both 5 and 8 GHz. 
Is this feature a counter- j et? In our images it is not as 
"jet- li ke" as in Figure 2 of Kollgaard, Ward le, fc Roberts 
(|1989|), an d it does not meet the criteria of |Bridle fc Per- 
ley |(|l984[ ). Nonetheless, the feature is real; at a distance 
ol 1 from the core the intensity is about 1 mJy/beam, a 

2 In those images that have been rotated in position angle by 
-50° we will refer to directions on the sky as they appear in the 
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Fig. 4. — Comparison of the ridge-line brightness in images made 
independently by two of the authors from the same 5 and 8 GHz 
data (Figure |2|). The black lines are at 5 GHz, the red lines are at 
8 GHz. Note the logarithmic intensity scale. 




Fig. 5. — Spectral index image of the quasar 3C 345 made from 
5 Sz 8 GHz data, shown over contours of total intensity at 5 GHz. 
The color represents spectral index on the scale to the top. The 
image was derived from Figure [2] and has been rotated by —50° 
in position angle. See Table [I] for image data and figure details. 

factor of about 20 fainter than the main jet at the same 
core distance, and so we take J = 20 as a lower limit on 
the jet-counter-jet ratio. In Figure [7] we show theoret- 
ical curves of constant jet-counter-jet ratio in the 
plane, assuming that the jet and counter-jet lie along 
the same line and taking a = —0.85. Examination of the 
curves shows that the fluid speed must be B > 0.50. Thus 



our images suppor t the conclusion of Kollgaard , Wardle, 
fc Roberts] ([1989) that at least mildly relativistic fluid 



speeds are required on kiloparsec scale s in 3C345. More 
stringent constraints are derived by Roberts & Wardle 
(2012) from a model of the polarization. We also show 
the curve j3 = cos i that expresses the constraint that the 
photons we observe are emitte d perpendicular to the je t 
axis in the frame of the fluid (Ro berts fc Wardle||2012 ). 



rotated images; note that the true position angles are 
-40°, "East" = +140°, "North" = +50°, and "South" 



'West" 
130°. 



3.2.3. Jet Opening Angle 

The shape of the jet may play an important role in un- 
derstanding the pec uliar electric vector configuration in 
this object (see £4.2 below). We examined the divergence 
of the main body ol the jet using a set of transverse slices 
taken from the image in Figure [3] We estimated the half- 
width at half-maximum of the jet in two ways, by fitting 
Gaussian profiles to the slices, and by searching for the 
location of the half-maxima on each side of the peak, 
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Fig. 6. — Intensity (5 GHz in black, 8 GHz in red) and spectral 
index (blue) along the ridge-line as a function of distance from 
the core, derived from the images in Figure pi The error bars 
correspond to an RMS uncertainty of 0.5 mJy/ 




Fig. 7. — Curves of constant jet-counter-jet ratio J in the (/3,i) 
plane. From left to right the ratios are 2, 5 10, 20, 50, 200, 2000, 
Sz 5000, computed for a spectral index of a = —0.85. The heavy 
line is the constraint that the photons emerge perpendicular to 
the jet axis in the frame of the fluid, (3 = cosi. 



and subtracting the beam quadratically. The results are 
shown in Figure [8j which shows deconvolved half-widths 
and fits. The linear least-squares fits yield half-opening 
angles of 9.31° and 9.43°, depending on whether or not 
the fit is forced to go through the origin. We adopt an 
apparent half-opening angle of <p a = 9.4°. 

3.3. Polarization Structure 
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Fig. 8. — Deconvolved half-widths and fits at 8 GHz; the data 
come from Figure [3] The data points are the means of Gaussian 
fits and simple searches for the half-power points, and the error 
bars are standard errors using these two estimates. The linear 
least-squares fits yield half-opening angles of about 9°. 

The linear polarizatiorjj structures at 5 and 8 GHz 
are shown in Figure |9{ a nd the fractional polarization is 
illustrated in Figure |10| We confirm the discovery by 

that the electric 
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vectors in the jet are not transverse to the jet but twisted 
at a significant angle Here we are able to resolve the 
polarization of jet in the transverse direction, and find 
that the twist is greatest, about 35°, down the center 
of the jet and vanishes at both edges. The fractional 
polarization of the jet is about 25% down the center of 
the main part of the jet, greater at the edges, and varies 
strongly with position in KB and HS. 

We first address the possibility that the twist in elec- 
tric vector position angles is due to Faraday rotation in 
the jet. This is ruled out as it is apparent from Fig- 
ure [9] that the electric vector distributions at the two 
frequencies are nearly identical. Polarization angle dif- 
ferences Ax = X5 — Xs across 3C 345 are shown in Fig- 
ure [TTJ where A% is displayed as position angle. An 
deviations from the vertical indicate Faraday rotation 
Everywhere in the source except the edges, |A%| < 10 
corresponding to about |RM| < 70 rad m -2 . The to- 
tal range of A% in the brighter parts, which is indepen- 
dent of absolute polarization position angle calibration, is 
max~^Xmin\ J$ 15°, or 100 rad m 2 . These amounts 
of Faraday rotation produce only small A\ at these fre- 
quencies, and we conclude that Faraday rotation is not 
the cause of the twist in the E vector directions in the 
main body of the je t. 



Rudnick & Jones (1983) found an integrated rotation 
measure of 29 ± 7 rad m -2 for 3C 345, but our measure- 
ment at the core yields —17 rad m -2 . Our absolute posi- 
tion angle calibration uncertainty of ±5° corresponds to 
a rotation measure uncertainty of about ±30 rad m -2 , 
marginally incompatible with the integrated measure- 
ment. We have examined data from the University of 
Michigan Radio Astronomy Observatory for this epoch, 
and they are compatible with our measurement of neg- 
ligible Faraday rotation. The source is highly variable 
and this is the likely cau se of the disagreement with the 
Rudnick fc Jones| ( |1983[ ) measurement. 

In highly-luminous radio sources with FR II structure 
such as 3C345 the inf erred magnetic field in the jets is 
typically longitudinal ( |Bridle et al. 1994), as evidenced 



3 The notation we use for linear polarization is as follows. The 
Stokes parameters Q and U are combined to make the complex lin- 
ear polarization P = Q + iU =pe 2ix , where m = (Q 2 + t/ 2 ) 1 / 2 // 
is the fractional polarization, p = ml is the linearly polarized in- 
tensity, and x — (1/2) arctan (U/Q) is the electric vector position 
angle. 

4 Between 4.86 and 8.44 GHz, A% = 14.6° for a rotation measure 
of RM = 100 rad m~ 2 . 
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Fig. 9. — Linear polarization images of the quasar 3C 345 with 
ticks showing the direction of the observed E vectors, made from 
(a) 5 & (b) 8 GHz VLA A-array data. The contours are of linearly 
polarized intensity. These images have been rotated by —50° in 
position angle. See Table ^ for image data and figure details. 




Fig. 10. — Color representation of fractional linear polarization 
over contours of total intensity at (top) 5 GHz and (bottom) 8 GHz. 
The color scales are 1000m. See Table [l] for image data and figure 
details. 

by orthogonal E vectors, while in low- luminosity FR I 
sources the inferred field is transverse (Bridle & Per ley 
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Fig. 11. — Linear polarization position angle differences in 3C 345 
shown over contours of linear polarized intensity at 5 GHz. Regions 
with no Faraday rotation will have vertical vectors. This image has 
been rotated by —50° in position angle. Data from Figure [9] See 
Table [l] for image data and figure details. 

1984|). Thus the linear polarization distribution in the 



straight portion of the 3C 345 jet is unusual in that the 
E vectors are neither longitudinal nor transverse to the 
jet axis. 

4. DISCUSSION 

4.1. Does the Jet in 3C345 Precess? 

The non-co-linear / structure of the jet seems to sug- 
gest that the jet was point ed in different directions at dif- 
ferent times in the past flKollgaard, Wardle, fc Roberts| 
1989). But the polarization structure strongly suggests 
that this is not the case, and instead that KB contains 
two impact points of a jet that is twice re-directed by in- 
teractions with the external medium, and that comes to 
an apparent end at the hotspot. We interpret the electric 
vector configurations at KBS, KBN, and HS as indicating 
compressed s heets of magnetic field seen roughly edge-on 
( |Laing||1980l ). In KBS and KBN the jet is redirected by 
interaction with the ambient medium, and at HS the jet 
ends in a hot-spot. This interpretation is consistent with 
a jet inclined at a small angle to the line of sight that is 
seen roughly from the side due to aberration because its 
fluid speed is highly relativistic. 

4.2. Origin of the "Twist" in the 3C 345 Electric 
Vectors 

Is the magnetic field in the 3C 345 helical, as it ap- 
pears to be? There are few similar examples in the lit- 
erature, the best observed one being M87. There the 
well-resolved jet shows a twisted filamentary structure 
in the radio (|Owen, Hardee, & Cornwell 1989) tha t has 



been ascribed to a Kelvin- Helmholtz instability (Hardee 
& Eilek 2011), and the magnetic field structure seems 
to follow the filament s. Optically it looks quite similar 
to th e radio image ( |Fraix-Burnet, Le Borgne, fc Nieto 
1989). In 3C345 we see no filamentary structure in the 
jet, and suggest instead that the twisted electric field pat- 
tern is the result of relativistic transfer effects. A helical 
magnetic field in a transparent homogeneous cylindrical 
jet would appear either transverse or longitudinal, not 
helical, due to cancellations between the back and the 
front (in a horizontal jet this means that Stokes U is 
identically zero). Two possibilities for producing a twist 
in the electric vector pattern suggest themselves; either 
the symmetry through the jet is broken by opacity or it 
is broken by relativistic effects. In the latter case, the 
gentle divergence in the profile of the jet demonstrated 
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in 1 3.2.3| suggests a similarly diverging velocity field. In 
such a velocity field, differential Doppler boosting due 
to differing line-of-sight velocity components in the front 
and back of the jet will break the symmetry, permitting 
there to be non-zero Z7, and thus creating a twist to the 
observed electric vectors if the underlying magnetic field 
is helical. It is shown in Roberts & Wardle (2012) that 
simple models for the magnetic field in 3C 345 coupled 
with differential Doppler boosting in a gently diverging 
jet can produce the "twist" seen in the electric vectors 
in the main body of the jet. This model can be used to 
place limits on the speed of the jet. 

5. CONCLUSIONS 

We have made sensitive VLA Array A-array images of 
the nearby bright quasar 3C 345 in total intensity and 
linear polarization at 5 and 8 GHz. The principle results 
of the analysis of these images are the following. 

1. 3C345 has a bright 7-8 Jy unresolved core with 
spectral index a ~ — 0.24 {I v oc v a ). 

2. The jet consists of a straight section of projected 
length about 2.3" (15 kpc), and two bright knots 
at the end of the jet that do not lie on the same 
line. 

3. The spectral index of the jet ranges over —1.1 < 
a < —0.5. The spectrum flattens significantly be- 
tween the end of the main body of the jet and the 
first knot at the end, and between that knot and 
the final one. 

4. The jet diverges slightly with an apparent half- 
opening angle of about 9°. De-projec ted, the jet 
opening angl e is significantly smaller ( |Roberts &; 



Wardle 



2012) 



5. There is a faint near-core feature whose base ex- 
tend from the core at (true) position angle ~ 55°. 
This feature does not appear to be a true counter- 
jet, but instead is part of an extended lobe seen in 
projection. The lack of a true counter-jet requires 
fluid speeds j3 > 0.5. 

6. The fractional polarization in the jet is high, rang- 
ing over 0.2 < m < 0.5, and is systematically 
greater at the edges. 



7. The non-co-linearity of the main body of the jet 
and the knots at its end seems to indicate that 
the jet was pointed in different directions at differ- 
ent times in the past. However, the polarization 
structure suggests instead that the knots represent 
three impact points with the external medium of a 
jet that is twice re-directed by interactions. 



3. The polarization behavior in the knots at the ap- 
parent end of the jet is consistent with the possi- 
bility that the jet, despite being at a small angle 
to the line of sight, is seen effectively from the side 
due to relativistic effects for a fast-moving flow. 



9. Unusually, the electric vector directions in the main 
body of the jet are neither longitudinal nor trans- 
verse, but are twisted. 
10. The "twist" in the electric vectors is not due to 
Faraday rotation, which is a few degrees or less 
everywhere, and is significantly less in the core, the 
brightest parts of the jet, and the knots at the end. 
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Image & Figure Properties for 3C345. 
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0.29 


0.36 


8.08 
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0.46 


0.30 


8.14 


1.0 
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0.40 


0.30 


7.11 


1.0 
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1.4 


0.20 


6.99 


1.4 
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0.46 


0.30 


8.14 


1.0 
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0.78 


0.30 


0.227 


0.71 
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0.99 


0.30 


0.209 


0.71 




1U 




0.78 


0.30 


0.227 


0.71 








0.99 


0.30 


0.209 


0.71 








0.78 


0.30 


0.227 


0.71 



NOTE: All contours are in steps of factors of a/2- 



